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Abstract 

Many nonlinear quantum phenomena of intense laser-atom physics can be intuitively explained 
with the concept of trajectory. In this paper, Bohmian mechanics (BM) is introduced to study a 
multiphoton process of atoms interacting with the intense laser field: above-threshold ionization 
(ATI). Quantum trajectory of an atomic electron in intense laser field is obtained from the Bohm- 
Newton equation first and then the energy of the photoelectron is gained from its trajectory. With 
energies of an ensemble of photoelectrons, we obtain the ATI spectrum which is consistent with 
the previous theoretical and experimental results. Comparing BM with the classical trajectory 
Monte-Carlo method, we conclude that quantum potential may play a key role to reproduce the 
spectrum of ATI. Our work may present a new approach to understanding quantum phenomena 
in intense laser-atom physics with the image of trajectory. 



* Electronic address: qycai@wipm.ac.cn 



I. INTRODUCTION 



In intense laser field (ILF), an atom may absorb multiple photons, more than that re- 
quired for ionization, and then it will eject a high energy photoelectron. Such a nonlinear 
phenomenon of multiphoton process is called above-threshold ionization (ATI) In 
general, the photoelectron spectrum of ATI consists of multiple peaks, separated by one 
photon energy . Many theoretical methods have been developed to study this multi- 
photon phenomenon, such as solving the time- dependent Schrodinger equation P, ?! , or the 
semiclassical trajectory methods including two-step model 



approach in the strong field approximation 



Bohmian mechanics (BM) 
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3, 



and Feynman's path-integral 



or called quantum trajectory method [15[, is an 



alternative to quantum mechanics. It has been successfully used to study some fundamental 



quantum phenomena such as tunneling [l6|, llTJ and the two-slit experiment [18|] . Recently, 
BM has been extensively applied to study many novel quantum processes in physics and 



chemistry, such as atom-surface physics 



dy man 



2l|, photodissociation of NOCl and NO2 



electron transport in mesoscopic systems 
and the chemical reactions {23]. It has also 



3een regarded as a resultful approach to studying chaos 
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26l . |27| and decoherence 



28|. 



In this paper, BM is introduced to intense laser-atom physics to study the multiphoton 
phenomenon of ATI. The trajectory of each atomic electron in an ensemble can be deter- 



lich is a subtle transformation of 
. Then the energy of an ionized 



ministically obtained from the Bohm-Newton equation w^ 
the Schrodinger equation and equal to it statistically [12 
electron, or called photoelectron, can be calculated from its trajectory. By rearranging all 
photoelectrons according to their energies, we gain a photoelectron spectrum of ATI which is 
consistent with the previous theoretical and experimental results. Finally, we briefly discuss 
the classical trajectory Monte-Carlo method (CTMCM) and conclude that a term in Bohm- 
Newton equation, called quantum potential, may play a key role to precisely reproduce the 
spectrum of ATI. 

This paper is organized as follow: We will briefly introduce BM first. Then we show the 
Hamiltonian for the hydrogen atom in ILF. The numerical solution of the time-dependent 
Schrodinger equation and details of BM for ATI together with the photoelectron spectrum 
are given. Next we compare our result with the previous theoretical and experimental results 
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of ATI. Finally, we compare BM with CTMCM and then conclude. 



II. THE QUANTUM TRAJECTORIES FORMALISM OF BOHMIAN MECHAN- 
ICS 

Bohm-Newton equation comes from a subtle transformation of the time-dependent 



Schrodinger equation [12|, ll3(]. A wave function ip can be expressed as ip{x,t) = 
t)e*'^^^'*)/'*, where R and 5* are real functions. Inserting ip{x,t) into the time- 
dependent Schrodinger equation, we obtain two equations by separating the time- dependent 
Schrodinger equation into real and imaginary parts. The real part gives 

|+(Z^ + v, + e.o (1) 

ot 2m 
and the imaginary part has the form 

§^ + V(p^) = 0, (2) 

where Q{x,t) = —^^^5 p{x-,t) = R^{x^t) and v = V S{x,t) /m. We find that equation 
(2) looks like the classical continuity equation, and equation (1) can be reduced to the 
classical Hamilton- Jacobi equation if the term Q{x, t) was ignored. In BM, Q{x, t) is usually 
called quantum potential and plays a crucial role for the appearance of quantum phenomena 
|l2 , ll^ . From the standpoint of classical mechanics, a Bohm-Newton equation of motion 
for a Bohmian particle can be constructed: 

md'^x/dt'^ = -V{V + Q). (3) 

In fact, according to the definition above, a much simpler equation of motion can be used 
to obtain quantum trajectory instead of equation ([3]): 

dx/dt = VS{x,t)/m. (4) 

Usually, we first solve the time-dependent Schrodinger equation to obtain ip{x, t), and hence 
S{x, t). Secondly we integrate equation (jlj) to obtain the quantum trajectory of the particle. 

Statistically, Bohm-Newton equation is completely equal to Schrodinger equation. All 
quantum phenomena described by Schrodinger equation can be reproduced by Bohm-Newton 
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equation, in principle, with a properly initial ensemble distribution. In particular, Bohm- 
Newton equation can help us to understand quantum world intuitively from the view of 
trajectory. In semiclassical method, the multiphoton phenomenon of ATI can be intuitively 
explained with the concept of electron trajectory. Therefore, quantum trajectory method is 
used to study the phenomenon of ATI in intense laser-atom physics in this paper. 



III. NUMERICAL SOLUTION OF THE TIME-DEPENDENT SCHRODINGER 
EQUATION 

In this paper, we study the system of hydrogen atom in ILF, using the one-dimensional 
model atom of hydrogen , I29I . I30I . The one-dimensional model, or namely the soft- 
core model, which has the long-range Coulomb tail characteristic of real atomic system, 
can approximate the system of hydrogen asymptotically. The Hamiltonian of the field-free 
atom is Hq{x) = — ^^^^^ s-nd the atom-laser interaction is Ht{x) = —xE{t), where 

E(t) is the laser field profile (atomic units are used throughout). Thus the time-dependent 
Schrodinger equation can be obtained i^^^^ = [Ho{x) + Ht{x)]4'ix,t). To solve the time- 
dependent Schrodinger equation, we first need to gain the eigenstates and eigenvalues of 
Hq{x): ipn{x) and En {n = 1,2, ...,N), respectively. Here we construct the eigenstates from 

M 

the B-spline basis, fn{x) = X] Cj'Bj{x) where Bj{x) is the B-spline basis and k is the 
order of B-sphne basis [3l|. The range of the variable x is confined to (— Xmax, a^max)- Using 
the standard diagonalization procedure, the eigenvalues En and values are obtained 
respectively. Secondly, the time-dependent wavefunction ip{x,t) is expressed in terms of the 
eigenstates as ^J^x^t) = ^ a„(t)v9„(x), where cu(t) is a time-dependent function. In this 
work we use Symplectic Algorithm method [32| to solve the time-dependent Schrodinger 
equation to obtain the coefficient a„(t), and hence the wavefunction ilj{x,t). 

Eo sin2(^) sin(cjt), < t < 3T 

Eq sin(co't), t > 3T 

where T = 2tt/uj, and Eq and u are the amplitude and angular frequency of the electric field 

in the laser pulse, respectively. Here Eq = 0.1 a.u. and u = 0.148 a.u. And we take: k = 7, 

M = 1200, = 1028 and x^ax = 600 a.u., with the time step of 0.0051824 a.u. The initial 

state iplx, 0) of the system is the ground state of the field-free one-dimensional model atom 

of hydrogen. 



In the present paper, the laser field profile is E{t) 
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FIG. 1: Quantum trajectories of electrons with different initial positions as functions of time. The 
initial positions xq are, from top to botton, 1.0 a.u., 0.9 a.u., 0.8 a.u., 0.7 a.u., 0.6 a.u., and 0.5 
a.u. 

IV. RESULTS 

After obtaining the time-dependent wavefunction ilj{x,t) and the companying S{x,t), we 
can numerically integrate equation (jlj) to get the quantum trajectory of electron with its 
initial position xq. Figure [T] shows the different electron trajectories with the corresponding 
initial positions as functions of time. Let's take the trajectory with the initial position 
Xo = 0.5 a.u. as an example: At the first stage before t = 450 a.u. about, the electron 
oscillates around the core, driven by the laser field. After that, the electron runs away from 
the core, following a linear line with small periodic oscillation. This implies that the electron 
has been ionized, i.e., a photoelectron has been ejected from the atom. We have obtained 
lots of trajectories of electrons with different initial positions Xq . All of them have the 
similar characters as explicitly shown in Figure [H We then calculate the average kinetic 
energy (Ek) of each ionized electron in one period of oscillation from its trajectory. Because 
of the long-tail potential of the model atom, the energy of the photoelectron is thus obtained 

E = {Ek) + V{x) 

atom 1 (5) 
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FIG. 2: The initial electron density distribution |'0(x,O)| in the range (—5,5). 

where V{x)atom is the long-tail potential of the model atom. In this way, we have shown 
the energy of a photoelectron can be obtained directly from the Bohm-Newton equation. 
In the traditional semiclassical trajectory methods, however, the energy of photoelectron is 
obtained from a semi-empirical formula The total energy of a free electron in the laser 
field is the sum of the ponderomotive energy and the translational energy: E = Ep + ^ {v)'^ , 
where Ep = {Eq/2uj)'^ and (v) is the electron translational velocity in the laser field. 

Next, we calculate lots of trajectories of atomic electrons in an ensemble to gain the 
energies of each photoelectron, and hence the photoelectron spectrum of ATI. Accordin g to 



Bohmian theory, the initial distribution of atomic electrons in an ensemble is {ipix, 0)p 12, 



ISj. In this paper, ?/'(x,0) is the ground state of the field-free one-dimensional model atom 
of hydrogen, which is obtained by numerically solving the time-independent Schrodinger 
equation: Hqi/j = Eip. Figure [2] explicitly shows the initial density distribution of the 
electrons. We have chosen 8,192 initial positions distributed in the range (—5, 5) with the 
distribution density 0)p (note that 99.9% of atomic electrons in the ensemble is in 
the range (—5, 5) at the time t = 0) and gotten the corresponding trajectories with the 
propagation time 16T. Then we gain the energies of the corresponding photoelectrons by 
equation ([5]). At last, by rearranging these photoelectrons according to their energies, we can 
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FIG. 3: ATI spectrum with three peaks at about lAiv, 2.4a;, and 3Auj, respectively. The inset 
is the ATI spectrum obtained by solving the time-dependent Schrodinger equation at the time 
t = 16.25r, after considering the ponderomotive after-acceleration of electron in intense laser field. 

gain a photoelectron spectrum of ATI, which is shown in figure [3] with three peaks at about 
1.4a;, 2.4a;, and 3.4a;, respectively. These peaks are contributed from different quantum 
trajectories with the corresponding initial electron positions Xq. Generally speaking, the 
photoelectrons with the initial positions near the core contribute to the first peak and the 
photoelectrons with the initial positions not near the core have the contribution to the 
second peak, the third peak, etc. For example, the photoelectron from the initial position 
Xq = 0.5 a.u. with the energy l.Sa; contributes to the first peak, and the one from xo = 1.8 
a.u. with the energy 2.4a; contributes to the second peak. 

Our result of ATI spectrum from BM is consistent with the previous experimental and 
theoretical results. In experiment, the positions of the peaks in a ATI spectrum are ac- 
curately given by the formula E = Eg + nu, where Eg is the ground state energy of the 
field- free atom and n is the number of photons absorbed by the atom P]. When Eg+n'uj > 
but Eg + n'uj — Ep < 0, due to the AC Stark shift, the channel for n'-photon absorption is 
closed and the corresponding peak in the ATI spectrum is suppressed, which is called chan- 
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nel closing and has been demonstrated experimentally 33|]. In our work, uj = 0.148 a.u. and 
the ground state energy of the model atom is Eg = —0.6697 a.u., so the peaks should appear 
at lASuj, 2.48co', and 3.48co', respectively. We have also reproduced the ATI spectrum by 
solving the time-dependent Schrodin ger equation [6]. After considering the ponderomotive 
after-acceleration of electron in ILF [3j|, we gain the ATI spectrum (the inset in figure [3]) 
which agrees with that from BM. In this way, we have shown that the ATI spectrum we 
obtained with BM is consistent with the previous experimental and theoretical results. 



V. DISCUSSION 



It is interesting that BM can well describe the process of ATI, in particular after compar- 
ing with the application of CTMCM to intense laser-atom physics. In the past, CTMCM has 
been used to study the process of ionization of the hydrogen atom in ILF 35|. It has been 
definitely pointed out that CTMCM itself cannot give the correct ionization rate, unless an 
additional tunneling process was considered (see figure 1 in 35|]). In CTMCM, motion of the 
electron is dominated by the classical Newton equation. In our calculation, however, Bohm- 
Newton equation has been used to gain the trajectory of the electron. We can thus conclude 
that the quantum potential Q{x,t) in equation ([3]) may play a crucial role to reproduce the 
spectrum of ATI in intense laser-atom physics. 



VI. SUMMARY 



In summary, we have used BM to describe the process of ATI in intense laser-atom 
physics. In our study, the energy of each photoelectron is obtained from Bohm-Newton 
equation and the obtained ATI spectrum is consistent with the previous theoretical and 
experimental results. Comparing BM with CTMCM, we conclude that quantum potential 
may play a key role to depict the quantum phenomenon. Therefore, our work may present 
a new approach to studying the quantum phenomena in intense laser-atom physics with the 
image of trajectory. 
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